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Abstract 

We analytically discuss the resonance conditions among several neutrinos in matter. The dis- 
criminant for the characteristic equation of the Hamiltonian is expressed by the coefficients of the 
equation. The result of the computation for the discriminants tells us that the neutrino energy 
and the matter density are in inverse proportion to each other at the resonance states in not only 



2- but also 3- and 4- neutrino models. 

m ■ PACS numbers: 14.60.Pq, 14.60.St, 13.15.+g, 96.40.Tv 

> ' 

o 

CO 

\o 

O 

o 

a,: 
6-;, 



X 



*Also at Radioisotope Center, Kyushu University, 3-1-1 Maidashi Higashi-ku, Fukuoka 812-8582, Japan; 
Electronic address: kamo@sci.kumamoto-u.ac.jp 



I. INTRODUCTION 

Neutrino oscillations mean transitions among neutrino flavors, and give an interpretation 
of neutrino phenomena l|, |2j, y, |4J found in the observations of the solar and the atmospheric 
neutrinos and the LSND experiment. The mass squared differences are some parameters of 
the neutrino oscillations, and are determined by some results of the observations and the 
experiments mentioned above. In order to explain the results of those within one framework, 
the three kinds of mass squared differences are needed. Existence of them means that four 
neutrinos have different masses from each other. It is a reason to discuss the 4-neutrino os- 
cillation. Four neutrinos are classified into active neutrino flavors (V e , v^ u T ), which interact 
with leptons in the weak interaction, and sterile neutrino (^ s ). The sterile neutrino does not 
have the weak interaction. 

Recent analyses of the experiments and the observations disfavor the four neutrino flavors. 
The possibility of an oscillation v e <-» v s is strongly excluded by the analyses. However, the 
result of the LSND experiment causes the maximum one of the three mass squared differ- 
ences, and gives grounds for the 4-neutrino models. The upcoming MiniBooNE experiments 
J5( may form a conclusion about this discrepancy. Whatever conclusion the experiment leads, 
it is useful to consider the neutrino oscillations with the sterile neutrino in the condition 
different from that of the experiment. 

The neutrino oscillation pattern in vacuum can get modified when the neutrinos pass 
through matter. This is known as the Mikheyev-Smirnov-Wolfenstein (MSW) effects[6], 
which can be described by an effective Hamiltonian. The Hamiltonian is expressed by 
the sum of the vacuum Hamiltonian and the interaction with the charged and the neutral 
— Q. In „ rder to write the 4 x 4 Harniltonran nratris in .natter, we introduce the 
4x4 mixing matrix which is an extension of the 3x3 Maki-Nakagawa-Sakata-Pontecorvo 
(MNSP) matrix|8|,|9|. 

Analytical calculations of active 3- neutrino oscillations in matter have been performed! 10. 

nfinn n 

L, 12. 13. 14J. In our previous papcqllj. we have given analytic expressions for the time- 
evolution operator in the 4-neutrino oscillation and the transition probabilities in the pres- 
ence of constant matter densities. 

In the results of our calculations, the value of the matter density which causes the neutrino 
resonance varies with the average of the neutrino energies. In general, the resonance occurs 



when two of the energy levels approach to each other. However, the calculation to get 
the energy levels is difficult because we have to solve the characteristic equation of the 
Hamiltonian. Then we propose to use the discriminant of the characteristic equation made 
from the Hamiltonian in matter in order to estimate the neutrino resonance conditions. In 
this calculation, we need not solve the characteristic equation. 

The outline of the article is as follows. In Sec. |Hl two kinds of bases to express neutrino 
states are introduced. These bases are connected to each other by a mixing matrix. To 
describe neutrino oscillations in matter, the effective Hamiltonian with charged and neutral 
currents is given. In Sec. IIII| we discuss the resonance conditions of neutrinos due to the 
characteristic equation of the effective Hamiltonian in matter. We then make use of a 
definition of the discriminant for an algebraic equation and a relation of the coefficients 
and solutions for the equations. In Sec. IIV1 the discriminant is concretely computed within 
a range expected that the neutrino mixings occur. And the graphical expressions of the 
discriminant are illustrated. Finally, we discuss the resonance conditions for the neutrino 
mixing in matter in Sec. 

II. FORMALISM 

We analyze the matter effects of 4-neutrino oscillation |Tj|, similar to those of 2- and 
3-neutrino oscillations. 

A. Two bases and a mixing matrix 

Neutrinos are produced in flavor eigenstates |z/ a )(a = e, /i, r, s). Between the source and 
the detector, the neutrinos evolve as mass eigenstates \u a )(a = 1,2,3,4). There are two 
kinds of eigenstates: \u a ) and \u a ). These eigenstates are defined by neutrino fields v a and 
v a corresponding to each eigenstate: z^|0) = \v), \u a ) = \a), \u a ) = \a), where a vacuum 
state is given by |0). In the present analysis, we will use the plane wave approximation of the 
fields. In this approximation, a neutrino flavor field u a is expressed by a linear combination 
of neutrino mass field v n : 
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v a = ^2U aa v a , (1) 

a=l 



"•<!• 



where U is a 4 x 4 unitary matrix with the elements U aa . If we write this relation in neutrino 
eigenstates, then 

l"> = !>:»• (2) 

An arbitrary neutrino state ^ is expressed in both the flavor and the mass bases as 



■0 



X) V«|a) = X^l a )' 



(3) 



a=e,n,T,s a=l 

where ^ Q an d if> a are the components of ip of the flavor eigenstate basis and the mass 
eigenstate basis, respectively They are related to each other in the form 



Ipa 



E UL1>a- 



(4) 



a=e,/i,r,s 



If we define the matrix elements as 
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(flavor) 



(« 
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(mass) 



u = (U aa ) 



M 



\ 



"02 



(5) 



(6) 



£4l U e 2 U e % U e 4 

Ufj,\ U^2 U/j,3 Upi 

U T \ U T 2 U T 3 Uri 

y U s i U s2 U s3 U S 4 j 

the relations (@J) of the components between the flavor and the mass eigenstates are rewritten 
in the matrix form, 

./(flavor) _ jj /(mass)^ /y\ 

The unitary matrix [/ is the mixing matrix of four neutrinos. There are 6 mixing angles and 
3 phases as parameters of U, in the case of 4 neutrinos. 
A parameterization for U is given by 



U — R34R24R14R23R13R12, 



(8) 



where R a t and R ab are the real and the complex mixing matrices, respectively [la]. The 
structure of these matrices is 

(Rabjij = Sij + (cos0 a6 - l)(S ai S aj + 5 bi 5 bj ) 
4 



+sin 9 ab (8 ai 8 bj - 5 bi 5 aj ) , (9) 

(Rob) y = $ij + {cos6 ab - l)(8 ai 8 aj + 8 bi 5 bj ) 

+sine ab (e- iA °»8 ai 8 bj - e iA <*8 u 8 aj ), (10) 

where the mixing of mass eigenstates \a) and \b) is described by 6 mixing angles 9 ab : #12, #13, 
#14, #23, #24, #34 and 3 phases A ab : A 13 , A 2 4, A 34 . In the case of leaving the CP violations 
out of consideration, i.e., A ab = 0, R ab = R ab , and the mixing matrix U is a real orthogonal 



matrix 



111- 



B. Hamiltonian in matter 

In the mass eigenstate basis, the Hamiltonian 7io participating in the propagation of 
neutrinos in vacuum is given by 

H = dis,g(E 1 ,E 2 ,E 3 ,E i ), (11) 

where E a (a = 1, 2, 3, 4) are the energies of the neutrino mass eigenstates \a) with mass m a : 



E a = ^m a 2 +p 2 . (12) 

Here and hereafter, we assume the momentum p to be the same for all mass eigenstates. If 

we assume |p| ^> m a , 

2 2 

771 III "^ a I I 1 ™ a /in\ 

E a ~ \p\ H — r- r ~ \p\ H 1 (13) 

where E is the average of the neutrino energies. 

The interactions between neutrinos and matter are described by two kinds of additional 
potentials. One is the interaction of the charged particles (electrons) and its neutrino u e : 

Vbc = V e = V2G F diag(iV e , 0, 0, 0). (14) 

The other is the interaction of the neutral particles(e.g., the neutron) and active neutrinos 

V NC = V2G F diag(-^iV n , -IjV n , -IjV n , 0), (15) 

where Gp, N e , and N n are the Fermi weak coupling constant, the electron number density 
and the neutral particle number density, respectively. Note that we assume the particle 



number densities to be constant throughout the matter where the neutrinos are propagating. 
The interaction term (J15J) can be separated into two parts as 

^nc = V n + V, (16) 

V n = V2G F diag(0, 0, 0, +^i\y , (17) 

V = V2G F diag(-^iV n , ~N n , -^N n , ~N n ). (18) 

The interaction terms V e and V n are added to vacuum Hamiltonian Hq for the propagation 
of neutrinos in matter. However, the vacuum Hamiltonian and the interaction terms are 
written by the mass and the flavor eigenstates, respectively. 

The mass eigenstate is convenient for describing the neutrinos propagation in matter. 
The Hamiltonian in the mass basis is 

?4r ss) = -Hq + U- x VccU + U- x Vk C U 

= H + U' l V e U + U- x V n XJ - A n I } (19) 

where I is the 4x4 unit matrix. The matter densities A e , A n and A are defined by 

A e = V2G F N e = A, (20) 

^, = -L GWV „ = i^. (21) 



However, the interaction terms U~ 1 V e U and U~ 1 V n U are complicated 11 ] . 

The flavor eigenstate is convenient for describing the interaction. The Hamiltonian in the 
flavor basis is 

^(flavor) = UHoU -l + Vcc + y NC 

= UHoU- 1 + V e + V n ~ A n I 

= UHoU- 1 + diag(A e , 0, 0, A n ) - AnI. (22) 

Although the term of a vacuum Hamiltonian is complicated, the interaction terms are simple 
for matter densities A e and A n . 

III. CALCULATIONS OF RESONANCE CONDITIONS 

Some analytical calculations show that the transition probabilities for the neutrino os- 



cillations in matter depend on the matter density 



0, 



llj |. In these analyses, the neutrino 



resonance occurs at the critical (matter) density. The critical density also changes with 
variation of neutrino energy. The lower neutrino energy reduces, the higher the density to 
observe the resonance of neutrino becomes. 

The resonance of neutrino occurs when the energy levels in the presence of matter ap- 
proach to each other. To discuss the resonance conditions, we consider the relations between 
the neutrino energy and the matter density. The neutrino energy is one of eigenvalues of the 
characteristic equation for the matter Hamiltonian. Therefore, the behavior of the eigenval- 
ues decides the neutrino resonance conditions. 

A. Discriminant of the characteristic equation 

To simplify discussion, we consider the Hamiltonian 7i^ avor ) in the flavor basis in (j22j) . in 
which the interaction terms are simple for matter density A. The degree of the characteristic 
equation of the Hamiltonian is determined by the number of neutrino flavors we suppose, 
e.g., the characteristic equation for 3-neutrino oscillations is a cubic one. The roots ai 
(i = 1, 2, 3, • • •) of the characteristic equation become more complicated as the degree of the 
equation increases. It is tedious to solve the characteristic equation of higher degree. 

In general, the discriminant Dt n \ of the characteristic equation 

x n + c n _ x x n - x + c n _ 2 x n ~ 2 + ■ ■ ■ + c = (23) 

with degree n, which has constant coefficients Cf. (k — 0, 1, • • • , n — 1), are expressed by the 
roots Oii with respect to x of the equation: 

D{n)= II {^-a 3 f. (24) 

l<i<j<n 

The discriminant Df n ) includes the difference between eigenvalues oii and (Xj, that is, a 
effective energy difference. It suggests that we can use the discriminant D( n ) to judge the 
neutrino resonance conditions. If D( n ) becomes small, we can conclude that the resonance 
occurs, because D( n ) depends on the square («j — <x,) of any energy differences. Therefore, 
some resonances may occur when one of («j — atj) has (local) minimum value. 

The discriminant D^ n \ is polynomials of the roots a,, of which degree is 2 x n C2, e.g., 
6 for the discriminant of the cubic equation. The discriminant D( n ) is also written as the 
symmetric polynomials for a«. The fundamental theorem of symmetric polynomials tells 



us that every symmetric polynomials in a, (i = 1,2,3, • • • ,n) can be represented by the 
following elementally symmetric polynomials: 

OLi + a 2 H \- a n = c n _i, 

a x a 2 + ai<y 3 H h « n -ia n = c n _ 2 , 

ai«2«3 H + a n - 2 a n -ia n = c n _ 3 , 

ai« 2 ■■■«„ = c . (25) 

The solutions on and the coefficients c^ for the characteristic equation are related to each 
other. Indeed, this relations between a, and Ck are given by (|25jl. In brief, the discriminant 
-D(n)) which includes the effective neutrino energy differences, can be represented by the 
coefficients of the characteristic equation for the Hamiltonian in matter without solving the 
characteristic equation. 

follows: 



For example, a discriminant Di 2 \ for a algebraic equation x 2 + Co = is calculated as 



D{2) = (oil - ce 2 ) = (a 1 + a 2 ) - Aa x a 2 

= -4c . (26) 

The concrete form of discriminants D/$< and Du\ are given as follows: 

D {3) = -27c 2 -W (forx 3 + Cl x + c = 0), (27) 

D (4 ) = 256c 3 - 27ci 4 + 144c ci 2 c 2 

-128c V - 4ci 2 c 2 3 + 16c c 2 4 

(for x A + c 2 x 2 + c x x + c Q = 0). (28) 

Note that, for the algebraic equation ()23|) with degree n, a term c n _ix n_1 with power (n— 1) 
of a variable x vanishes by a suitable transformation: x — > x — c n _i/n. 

B. Traceless Matrix Tjy 

In order to find the explicit form of the characteristic equation, the Hamiltonian in the 
matrix form is separated into the diagonal and the traceless matrices. An arbitrary N x N 

8 



matrix M can always be written as 

M = T N + ^(trM)I N , (29) 

where TV and In are N x N traceless and unit matrices, respectively. Note that trTjv = 0. 
The characteristic equation for the matrix T/v is written by an eigenvalue A of T/v as 

= det(T N - XI N ) 

= \ N + cjv^A^" 1 + • • • + c 2 A 2 + ciA + c , (30) 

where the coefficient c is the determinant of T N , e.g., c = detT N and c^ (k — 1,2- ■ ■ N—l) 
are expressed by the sum of cofactors of the diagonal elements of T/v. But cn-i = 0, because 
Cjv-i is equal to the trace of the traceless matrix T N . 

In the 4-neutrino model, the Hamiltonian in the matrix form is given in (J19)) and (|22|) . 
Ignoring the CP violations, the trace of the matrix 7i^ lavor ' ) in ()22|) is 

tr ^ avor ) =E 1 + E 2 + E 3 + E, + A e - 3A„, (31) 

where we use the unitary conditions, e.g., U e \ '+U e 2 2 +U e3 2 '+U e i 2 = 1. Then, the Hamiltonian 
in matter is described as 

^(flavor) = ^ + 1 ( tr ^(flavor) j ^ (32) 

The characteristic equation for matrix T is written as 

A 4 + c 3 A 3 + c 2 A 2 + Cl A + c = 0. (33) 

The coefficients q (/ = 1, 2, 3, 4) are expressed in the matrix form as 

Co = detT, 

C\ = — C0fT(i) - cofT( 2 ) - C0fT( 3 ) - CofT( 4 ), 
C 2 = COfT(i2) + COfT(is) + COfT(H) 

+ C0fT {23) + C0fT (24 ) + C0fT(34), 

c 3 = -trT = 0, 

where the cofactors cofT( p ) (p = 1,2,3,4) of diagonal components T pp of T 4 and cofT( rs ) 
(1 < r < s < 4) of T rr and T ss are determinants of 3 x 3 and 2x2 matrices, respectively, 



e.g. 



cofT/ 



(2) 



E 



~PlP2P3P4 



Tip 1 d2 P2 T 3p3 T i 



4p 45 



coiTi 



(13) 



Pl,P2,P3,P4=l 

4 

/ , 6 PlP2P3P4 

P1,P2,P3,P4=1 



"lpi ^2p 2 "3p 3 ^4p 4 • 



The coefficients Ck of the characteristic equation for the matrix T4 are obtained from the 
components T rs of T 4 , and the discriminant _D( 4 ) is derived by these coefficients from (|2%|). 
The discriminants of T 2 and T 3 in 2- and 3-neutrino models, respectively, are given in the 
similar method as in 4-neutrino model. 

C. Resonance conditions for 2-neutrino oscillation 



The discriminant of the characteristic equation for the matter Hamiltonian is derived 
by the calculations shown below. For the simple case, a calculation for the discriminant of 
2-neutrino Hamiltonian is performed. 

The matter Hamiltonian for v a and vp is described as 



^(flavor) = un {va a ss) u -l + y^ + ^ 



NC 



E-^cos26 ^f- sin 26 



^ sin 29 



E + ^cos 26 



+ 
where 9 is the mixing angle, and 

E21 = -E-2 



A S 

/1 c u ea 




A n 5 sj 3 



(34) 



Ei, E 



5<* 



Eo 



(35) 



The interaction term in (}Mj) includes the factors A e 5 ea and A n 5 s p. These factors appear if 
we assume the electron neutrino v e and the sterile neutrino z/ s , respectively. In the case of 
2-neutrino Hamiltonian for v e and v„, we can leave A„5„a out of consideration. 



'M> 



n<J s /3 



From (|2*H|) . the Hamiltonian (pM|) is rewritten for the 2x2 traceless matrix T 2 : 

^(flavor) = T2+ l (tr ^(flavor) )/2) 
trH (flavor) = 2E + ^^^ + ^ 



(36) 

(37) 
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The concrete form of the traceless matrix T 2 is 



To 



E21 
~2~ 



cos 26 sin 26 
sin 26 cos 26 



:{A e S e 



j4 n <y fl/ ?)diag(l, -I) 



(38) 



The characteristic equation for the matrix T 2 is expressed as 

= det(T 2 - \h) 

= A 2 - - [(A e 5 ea - A n 5 s p - E21 cos2#) 2 + E 21 2 sin 2 26 

Then, the discriminant Di 2 ) is presented as follows 

D (2) = (A e 5 ea - A n 5 sP - E 21 cos 26) 2 + E\\ sin 2 26. 



(39) 



Therefore, the condition for the neutrino resonance, where f|39|) gets (local) minimum value, 
for two neutrinos u a and up is 



A e 5 ea - A n 5 s/3 = E 21 cos 26. 



Note that energy differences E 2 % are approximately given by 

E 2 \ ^ ej 2 — fj\ 



m 2 2 — mi 2 _ Am 21 



(40) 



(41) 



2E 2E ' 

where (|13|l has been used, and Am 21 is a squared mass difference. 

From (|4(Jj) and (|41|) . we conclude that, for the 2-neutrino resonance conditions, the average 
E of the neutrino energies and the factor (A e 5 ea — A n 5 s p) of the matter densities are in inverse 
proportion to each other: 



E(A e 5 ea - AJ, 



sj3) 



Ami 

2 ±cos26. 



(42) 



Note that (jj2j) gives a usual expression jig for the resonance condition if we set A n 5 s p = 0. 

IV. NUMERICAL ANALYSES OF RESONANCE CONDITIONS 

We apply some results obtained in the previous section to the 2-, 3- and 4-neutrino models. 
A discriminant explicitly contains the neutrino energy differences E a \, = E a — Eb (a, b = 
1,2,3,4 in the 4-neutrino model) in vacuum. These energy differences are approximately 



11 



given by the matter mass differences Am^ = m 2 — m 2 -, which are well-known quantities 
in the various neutrino oscillation experiments. In the similar way as showed in ([41^. 

Am ab 2 



E 



ab 



(43) 



2E ' 

where E is a average of the neutrino energies, i.e., in four- neutrino model, E = (E\ + E2 + 
E 3 + £ 4 )/4. 

These days, three kinds of mass squared differences are known. They are represented by 
Am 2 olar , Aml tm and A?til SND , which are used as the parameters in the solar and atmospheric 
oscillations and the LSND experiment. Using these mass squared differences, several distinct 



a 



types of mass patterns are possiblel?]. The phenomenology and the mixing matrix depend on 
the type of the mass scheme. We concentrate on the discussion on the (3 + l)i - scheme \ 111]. 
which is one of the several mass patterns in Fig. ^ This mass-scheme includes the usual 
3-neutrino mass pattern. 

In the mass pattern (3 + l)i, there are three close masses and one distinct mass. Let m 4 
and Am| 3 be the distinct mass and the largest mass squared difference, respectively. Three 
kinds of the neutrino mass squared differences are put as follows [17[ : 



Am 2 21 
Am 2 32 



Amf olar ^7xlO"°eV^, 



,2 

'atm 



Am: 



5 p ^ -2 
3^r2 



Am 2 



41 



Am 



LSND 



3 x l(T J eV 
-leV 2 . 



(44) 
(45) 
(46) 



Then, the energy differences E a b are expressed by 



E. 



21 



IE 



E: 



:-!2 



E 



Am\ x 



11 



IE ' 
E42 — E41 — E 2 i, 



E31 ■ 
E43 



E32 

■ E 41 



AW3 2 

IE ' 
1 E 2 i, 
- E 3 i, 



(47) 



where we suppose that the average E of the neutrino energies is 10 GeV or 10 MeV. 

n 

Next, we consider the approximate mixing matrix for the (3 + l)i-scheme |18i |: 
4=cose cos 5 4jCose cos<5 cose sin 5 sine ' 



h sin 5 



_i 

2 



72 cos6 ° 
^sm5 -j^cosS 



i-|sin5 

_i 
2 



V V2 



1 sine cos 5 ^isinecos<5 —sin S sine cose 



v^ k 



(48) 
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where 5 = S^n and e are small: < e < 0.1. Here the 3x3 sub-matrix that describes the 

loU 



mixing of the three active neutrinos has the bi-maximal form. The mixing matrix of ()48j) is 
given from (JBJ), © and (fTU|) by taking 

#12 = J, #23 = J, #13 = ^ 
#14 = 6, #24 = 0, 6*34 = 0, 

Ai 3 = A 14 = A 24 = 0. (49) 

In 2- and 3-neutrino models, we can also discuss the resonance condition by a sub-matrix 
of 4 x 4 mixing matrix. For example, in order to consider the 3-neutrino model, we may 
ignore the fourth elements in (J48J1 . taking e = 0. In this case, 4x4 mixing matrix (}4*%|) is 
handled as 3 x 3 matrix and a unity element. 

In this setting, we can get a numerical expression of the discriminant D^) and -D(4) as a 
smooth function of A: 

L>( 3) ~ -= [6.5 x 1(T 21 - 8.4 x 10~ 18 (AE) 
E 

+5.2 x 10- 12 (AE) 2 - 6.8 x W~ 9 (AE) 3 
+2.3 x 1(T 6 (AE) 4 ], (50) 

£> (4 ) ~ -=& [4.2 x 10~ 22 - 1.3 x 10~ 19 (AE) 



E 



+8.1 x 10~ 14 (AE) 2 - 1.0 x 10~ 10 (AE) 3 

+3.5 x 10~ 8 (AE) 4 + 7.1 x 10- 8 (AE) 5 

-3.2 x 10" 8 (AE) 6 - 1.4 x 10" 7 (AE) 7 

-3.4 x 10- 8 (AE) 8 + 7.1 x 10- 8 (AE) 9 

+3.5 x 10~ 8 (AE) W 

+1.7 x 10~ 18 (AE) n ], (51) 



where e = 0.1. The derivative —rP^ of D( n ) may be useful to estimate the neutrino resonance 
conditions. The solutions of the equation ,^" ) = will give the the turning value of the 



discriminant D( n ) ■ Although, analyses of the derivative of the discriminant _D( n ) for 3- and 



dD 



4-neutrino models are complicated, it is expected that the solution of the equation J?' = 0, 
that is a resonance condition, is given as an equation for A ,E, and the function of some 
parameters: 

AE = f(Am 2 ab ,9 ab ,A ab ). (52) 

13 



If the value in the right hand side of the (|52|) is fixed, e.g, as in (JUJ, P°T) - P^ an d dlSl) ; we 
can conclude that £? and A are in inverse proportion to each other. 

A. Normalized discriminant 

In this analysis, the discriminant of the characteristic equation for the Hamiltonian in 
matter has a small value. The value of the discriminant in vacuum is given by setting the 
matter density A to 0. Here we assume that the electron number density N e is equal to the 
neutral particle number density N n : N e = N n which means A e = A and A n = ^A in (|2(J|) 
and (EH). 

The effective energy differences |A a — Xb\ (a,b = 1,2, 3, 4, a ^ b), which are the dif- 
ferences between eigenvalues A a of the Hamiltonian in matter, are functions of the matter 
density A. If we set A = 0, then the effective energies A a can be regard as the eigenvalues 
E a of the Hamiltonian in vacuum: 

\X a - X b \ = \E a - E b \ ioiA = 0. (53) 



In a 2-neutrino (z/ e , v^) model, from (|41j) . we get a rough estimation for the discriminant 



D(2): 



D (2)\ A =o = ( A 2 - Ai) \ A=0 = (E 2 -E 



2| _ /-ei z^ \2 _ 



Am 2 i 21 " 



Am 2 



solar 



2E 



IE 

l(r- :, [<>Y 2 ] (for E = 10 lu eV). 

(54) 



2 

29UA/21 /r_ 7-1 _ 1n l0 



It may be easy to expect that, in the 3- and 4-neutrino models, the discriminant -D(3), -D(4) 
in vacuum have very small value, since the discriminant are defined by the product (J2IJ). 
Indeed, D {3) ~ l(T 82 [eV 6 ] and D {4) ~ l(r 142 [eV 12 ] for A e = A n = 0, e.g., in vacuum, by our 
rough calculations. Here we set E = 10 GeV. 

In order to concentrate to the argument about the behavior for D^ as the function of 
A, we define the normalized discriminant D n which is expressed as the value of the D n (A) 
divided by that for A = 0: 
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Hereafter, D4 means the normalized discriminant for the 4x4 matrix which is the operator 
of the 4-neutrino (v e , z/ M , u T , v s ) model. Moreover, D3 specializes the normalized discriminant 
for the 3x3 matrix of Hamiltonian, in which we deal with the three active neutrinos: 
(y e , Up, v T ). In order to compare the 3- and 4-neutrino models with the usual neutrino 
phenomenology, we consider the two neutrinos (z/ e , v^) for D 2 . 

B. Graphical expression of resonance conditions 

As an illustration of the resonance phenomena, the normalized discriminants D n (n = 
2,3,4) are plotted as a function of the matter density A in Fig. |21 Here we set the average 
E of the neutrino energies to 10 GeV, which is a value of neutrino energy used in previous 
analyses |l(llll|. 

From Fig. |21 the following results can be derived. If the mixing angles #13 and #14 are 
small, the third- and fourth-neutrino effects appear as a resonance. The sharp drop in the 
Fig. 12 is a good reflection of the neutrino resonance. Indeed, the surviving probabilities P ee , 
which are discussed in jlOj and [llj, for neutrino oscillations in matter have sharp drops. 
The smaller a mixing angle reduces, the sharper the drops for the discriminant D( n ) and the 
surviving probabilities become. 

One effect for the third neutrino v T appears beyond A ~ 10~ 13 eV. The other effect for 
the fourth neutrino v s also occurs beyond A ~ 10 _10 eV. It means that the matter density 
to cause the resonance for v e <-► v s or v e <-> v T varies with the squared mass differences: 

1(r '" l ^ Am LSN D | (56) 



10- 13 Am!? 



aim 



The above results are argued for E = 10 GeV. The normalized discriminant for E = 
lOMeV is shown in Fig. El The lines in Fig. [2] and Fig. El have the similar shapes except 
for the value of the matter density A. The difference of two results occurs due to the 
difference of the average E of neutrino energies. For instance, one resonance for v e <-» v T in 
Fig. |21 (E — 10 GeV) occurs for A ~ 10~ 13 eV, the other resonance for v e <-> v T in Fig. El 
(E = lOMeV) occurs at A ~ 10 _10 eV. Therefore the average E of the neutrino energies 
and the matter density A are in inverse proportion to each other as shown in (J4*2"|) . i.e., 



EA = =^tsi cos 2 13 ~ 10~ 3 [eV 2 ]. (57) 
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Note that we set the matter density to A e = A, A n = (1/2) A In the case of the resonance 
v e <-► u s , 

EA = Am 2 LSND cos26 u ~ l[eV 2 ]. (58) 



V. DISCUSSION 

The main result of our analysis is derivation of the discriminant of the characteristic 
equation for the Hamiltonian in matter, which is discussed in flavor basis. The concrete 
form of the discriminant for 2-, 3-, and 4-neutrino models is given in ([26)1. (|27|). and (|28|h 
respectively. The effective Hamiltonian in the 2-neutrino model is also calculated concretely. 
The discriminant made from the effective Hamiltonian is a function of the matter density A 
and the average E of the neutrino energies. The calculations for 3- and 4-neutrino models 
are straightforward but the forms of the discriminants are complicated. One of the numerical 
expression of the discriminant is given in (joTIj) and (|51|). The graphical expression of the 
discriminant are also shown in FigJ2]and Fig|3J 

One of the resonance conditions for the neutrino oscillation is described in ()42|) . It shows 
that E and A are in inverse proportion to each other, as confirmed in Fig. and Fig. 
01 That is, the lower the neutrino energy reduces the higher the density to observe the 
resonance between the neutrinos becomes. The concrete form of the resonance condition for 
2-neutrino model is shown in ()57|) and (|58|). From (|5Uj) and 1)51)1. the discriminant _D( n ) is a 
smooth function of A The derivative A n) of D< n ) m ay be useful to estimate the neutrino 



resonance conditions. Although, analyses of the derivative of the discriminant Di n ) for 3- 
and 4-neutrino models are complicated, we can expect that the neutrino resonance condition 
will be given as ()52j) . 

The conditions that the discriminant Di n ) has a local minimum value give the resonance 
conditions for neutrino oscillation. In our analysis, it is expected that the resonance of 
v e *^ v T and v e <-► v s occur at A = 1CT 13 eV and A = 10 -10 eV, respectively, for E = 10 GeV. 
This results can be confirmed by the previous calculations [lOj, |jj| in 2-, 3- and 4-neutrino 
models, respectively. For instance, the surviving probabilities P ee for 4-neutrino oscillation 

is given as 

34 . [T. A™ , 2 1 

(59) 



i" 4 E J2 B a B b sm 

a=l b=a+l 



L Am ab 21 



2 IE 
where B a is a factor determined by the mixing matrix and the characteristic equation of the 
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matter Hamiltonian in mass basis 111 ]. The factor L is a baseline length which stands for the 
neutrino flying distance. Fig. 0]is a result of the recalculation of the transition probabilities 
for new conditions which we have used in this paper. One can find that the value of the 
matter density A that the resonance occurs in Fig. |2] coincides with the value of A that the 
sharp drop of the transition probabilities occurs in Fig0] 

Is the condition to observe the neutrino resonance realistic? The electron matter density 
A is 1 x 10~ n eV at the center of the sun. The average of the electron matter density in 



earth is about 5 x 10 13 eV 



0. 



In the case of E = 10 GeV, it may be possible to observe 
the resonance between v e and v T since the condition is A ~ 10~ 13 eV which corresponds to 
the density of the earth. But in other analysis it may be difficult to observe the resonance. 
One may consider that the resonance will occur in the high energy and low matter density, 
which is restricted by (}4*2*j) . But there is a limit of the neutrino energy to observe the 
neutrino oscillation. From ()59|) . the baseline length L to observe the neutrino oscillation 
and the average E of the neutrino energies are proportional to each other. For instance, 
the condition to observe a neutrino transition v e — > v s for E = 10 GeV and E = 10 MeV is 
L > 25000 m and L > 25 m, respectively, where we assume the condition 

L = Am ob 2 -L 

A/4 2nE K ' 

Then, we have considered that the neutrino oscillation can be observed when the neutrino 
flying distance is more than A/4, which A corresponds to the wavelength in (JSTJj) . 

The 3-neutrino model with a small and bi-maximal mixings and the 4-neutrino model 
will be tested more closely by the upcoming experiments. In this paper, the resonance 
conditions for 2-, 3-, and 4-neutrino mixings are given without a preconception. Our result 
of the resonance conditions, which are given by using the discriminant, does not deny the 
possibility of the 4-neutrino oscillation, but the observation of the resonance is not realistic. 
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FIG. 1: Several mass patterns for 4-neutrino schemes. 
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FIG. 2: Normalized discriminants D n (n = 2,3,4) as a function of the matter density A for a 
(3 + l)i-scheme, where 9 12 = f , 23 = f , #13 = ifo 71 "' #14 = ifo 71 "* #24 = 0, #34 = 0, A i3 = A 24 = 
A 34 = 0, Amli ~ 7 x 10" 5 eV 2 , Am§ 2 ~3x 10" 3 eV 2 , Ara^ ~ leV 2 , E = 10 GeV. The solid, 



broken and dotted lines are normalized discriminants D4, D3 and D 2 , respectively. 
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FIG. 3: Normalized discriminants D n (n = 2,3,4) as a function of the matter density A for a 
(3 + l)i-scheme, where 9 12 = f , 23 = f , #13 = ifo 71 "' #14 = jfo 71 ": #24 = 0, #34 = 0, A i3 = A 24 = 
A 34 = 0, Amli ~ 7 x 10~ 5 eV 2 , Am| 2 ~ 3 x l(T 3 eV 2 , Am| ~ leV 2 , E = 10 MeV. The solid, 



broken and dotted lines are normalized discriminants -D4, D3 and -D2; respectively. 
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FIG. 4: The surviving probability P ee of the electron neutrino transition as a function of the matter 
density A for the (3 + l)i-scheme, where 6*12 = f , 6*23 = \ , #13 = jfo 71 ") $14 = e > #24 = 0, #34 = 0, 
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A 34 = 0, Am^ ~ 7 x 10~ & eV 2 , Am| 2 ~ 3 x lO^eV 2 , Am^ ~ leV 2 , £ = lOGeV 



and L = 10000 km. The solid and broken lines show the transition probabilities for e = j|q7t,0, 
respectively, which correspond to the 4-neutrino and the 3-neutrino cases. 
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